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CN . Abstract. 

'^ i The raise and peel model of a one-dimensional fluctuating interface (model A) is extended by 

1-^ I considering one source (model B) or two sources (model C) at the boundaries. The Hamiltonians 

'sj" ' describing the three processes have, in the thermodynamic limit, spectra given by conformal 

^^ ■ field theory. The probability of the different configurations in the stationary states of the three 

CN ■ models are not only related but have interesting combinatorial properties. We show that by 

.^ ! extending Pascal's triangle (which gives solutions to linear relations in terms of integer numbers), 

Psj I to an hexagon, one obtains integer solutions of bilinear relations. These solutions give not only 

O ■ the weights of the various configurations in the three models but also give an insight to the 

(^ ■ connections between the probability distributions in the stationary states of the three models. 

■^ . Interestingly enough, Pascal's hexagon also gives solutions to a Hirota's difference equation. 

§: 

^; 1 Introduction. 

c^ ' Recently much interest have been devoted to surprising appearence of the ASM (alternating 

C . sign matrices) combinatorics in the properties of the ground state wave function of XXZ 

>■ ! spin chain at a particular value A = —1/2 of its anisotropy [1]. In subsequent investigations 

'k>( I [2]-[12] a number of models having ground states with interesting combinatorial properties 

^ I was found, including the dense 0(1) loop model, the rotor model and the "raise and peel" 

model of fluctuating interface. Two important facts are common for these models. First, all 

they admit a purely algebraic description in terms of an appropriate version of a (quotient 

of) Temperley-Lieb (TL) algebra (see [10]). Secondly, the TL algebras always come in a very 

specific semigroup regime. This latter fact was used in [5, 10] to interpret the loop model 

as a stochastic process and thus, to give a physical interpretation to the components of the 

ground state wave function as (unnormalized) probabilities of various configurations. 

In the present paper we extend the raise and peel model^ (hereafter called model A) 
obtained for the TL algebra to two other cases. In order to do so, we consider the boundary 

^E-mail: pyatov@thsunl.jinr.ru 

^For detailed discussion of physical properties of the raise and peel model the reader is referred to papers 
[12, 13]. 
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extension of the TL algebra which is called the blob algebra [14, 15, 11] (see Section 2). In this 
way, in the stochastic model one can introduce a source at one boundary (model B) or sources 
at the two boundaries (model C). The two new models are described in detail in Section 3. 
As shown in Ref.[16], in the continuum limit, the spectra of the two Hamiltonians giving the 
time evolution of the stochastic processes are given by characters of N=2 superconformal 
field theory. Here we are going to consider the combinatorial properties of the stationary 
states only. 

In Section 4, based on numerical studies of small chains we make a series of conjectures 
for the weights of various configurations observed in the stationary states of the three models. 
In Appendix A we define Pascal's hexagon. We think that the content of this Appendix is 
interesting on its own. Using Pascal's hexagon one can get in a simple way the numerical 
results obtained for the stationary states for finite systems (see Section 4). 

The Pascal's hexagon is connected in a profound (although not yet understood way) to 
alternating sign matrices with various symmetries and as discussed in Appendix A with the 
solutions of a discrete Hirota's equation. 

After this work was almost completed, we learned from Jan de Gier that part of our 
results were obtained independently by Mitra et al [17]. We have also verified, as discussed 
in the Appendix, that Pascal's hexagon can be used to obtain properties of the stationary 
states of stochastic processes related to the periodic TL algebra discussed in [17]. 

2 Boundary extended Temperley-Lieb algebra 

2.1 Definition. 

We start with the type A Temperley-Lieb (TL) algebra [18] which as it is well known (see 
[19]) stands behind the Uq{sl2) symmetric XXZ spin chain. For the chain of L particles 
the corresponding TL algebra is generated by the unity and a set of (L — 1) elements e,, 
i = 1, . . . (L — 1), subject to relations 



e: 



2 



(g + g-i)e, , (2.1) 



CiCj = CjCi , V i, j : |i - jl > 1 . (2.2) 

Here parameter q G C\{0} of the algebra is related to the anisotropy parameter A = — (g + 
q~^)/2 of the spin chain. 

A boundary extension of this algebra is achieved by adding two more generators /o and 
fi together with the relations 

ei/oci = ei , ejfo = foej, Vj > 1, (2.3) 

eL-ifiei-i = ei-i , ejfi = fiej , V j < L - 1. (2.4) 

/o=«/o, fl = afL, fofL = fLfo, (2.5) 



where a, a G C. The algebra with one boundary element (either /o, or fi) called the blob 
algebra was analyzed in [14, 15]. The extension of the TL algebra with the two boundary 
generators was introduced in [11]. Unlike TL and blob algebras it is infinite dimensional 
and we are going to further extract its finite dimensional quotient. To this end we consider 
a pair of unnormalized projectors X^ and Y^. They are defined differently depending on a 
parity of L 



for L even: X^ 



L/2-1 
k=0 



Yr :-- 



(L-l)/2 

for L odd: Xl := /o H ^2fc , >l - 

fc=i 

In terms of these projectors reduction conditions read 



L/2-1 

/o 11 e2fc /l ; 

fc=l 

(L-3)/2 

n e2fc+l /l 
fc=0 



(2.6) 



X^YiXi = 6 Xi , YlXlYl = b Yl 



(2.7) 



The resulting quotient algebra is that one we shall further call the boundary extended TL 
algebra. It is finite dimensional^ and it depends on four parameters g, a, a and b E C. 



2.2 Graphical presentation. 

There are at least two ways in which the boundary extended TL algebra can be visualized. 
First one is a straightforward generalization of the diagrammatic realization of the blob 
algebra presented in [14]. For this one uses familiar "lines and loops" diagrams for the TL 
generators Cj and realizes the boundary generators /o and /^ as two different blobs lying, 
respectively, on the leftmost and the rightmost lines of the diagram (see [14]). Equivalently, 
one can draw boundary generator /o (/l) as a half-loop connecting the leftmost (rightmost) 
line to the boundary (see [11, 17]). 

The second way which we are using throughout this paper is the one suitable for modelling 
of growing interfaces (see [12]). One draws the TL generator Cj as a tile whose diagonal 
is lying on a vertical line with coordinate i and whose left and right vertices are placed, 
respectively, on vertical lines with coordinates i — 1 and i + 1. The boundary generators /o 
and fi are drawn as half-tiles with their longest sides lying on vertical lines with coordinates 
and L, respectively (see Figure below). 
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^Dimensions of the boundary extended TL algebras are calculated in an Appendix to [20] . 



The (half-) tiles can freely move along vertical axes unless they meet their neighbors. As- 
suming attraction forces acting among the (half-)tiles one represents word in the algebra as 
a collection of dense polygons built from the (half-)tiles and satisfying following conditions. 
All polygons are placed between vertical lines with coordinates and L and no vertical line 
lying between these two boundary verticals crosses the borders of (one or several) polygons 
in more then two points. 



2.3 The ideal X^. 

Of our main interest is the left ideal in the boundary extended TL algebra generated by 
Xl^. We denote this ideal as Xl. 

Consider graphical realization of a typical word in the ideal. As we have different defini- 
tions of Xl depending on a parity of L, separate pictures for the cases of L even and L odd 
are given below. 



L even: 




L odd: 




(2.9) 



Here we adopt a convention that multiplication from the left by elements Cj, /o, or fi amounts 
graphically to dropping their respective (half-)tiles up-down. Components of Xl are shown 
hatched on the pictures. 

As it is obvious from pictures (2.9) each word w in the ideal Xl is uniquely defined by a 
shape of the upper border h{w\x), < x < L, of its corresponding polygon (drawn in bold 
lines on the pictures). In turn, the border line is suitably encoded by its values at the integer 
points h{w\i) := hi{w). i = 0, 1, . . . , L. Assuming the height of the tile (= the length of its 
diagonal) equals 2 and taking the middle line of the bottom row of tiles as a reference axe 
one gets following prescriptions for a set of {/ii}i=o,i,...,L 



a) hi+i — hi = ±1 and hi > , V i ; 

b) hL is an even integer; 

c) there exists i such that hi G {0, 1} . 



(2.10) 



Here prescription b) results from our choice of Xl as an ideal generating element. With the 
choice of Yl one would constrain /ig to be even. Prescription c) arises from the reduction 
conditions (2.7). 

^Note that in case & 7^ the element Yl generates an isomorphic ideal. 
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In Ref.[21] a set of data {/ij}i=o,i,...,L satisfying conditions (2.10) is named an Anchored 
Cross path, L is called a length of the path. Anchored Cross paths of length L label effectively 
words in the ideal X^. There are 2^ different Anchored Cross paths of length L (for the proof 
c.f. Appendix of Ref.[20]) and thus, dimX^ = 2^. It is remarkable that the dimension of Xj^ 
coincides with the number of states of the chain of L spin=l/2 particles. This is not just 
a coincidence and the ideal X^ can be used for representation of an open XXZ chain of L 
spin=l/2 particles (see [16]). 

In considerations below we will use besides the set of Anchored Cross paths a pair of its 
subsets (or, equivalently, two subspaces in the ideal Xx,). Their definitions are given below. 

Ballot paths are the paths (2.10) with fixed endpoint Hl = 0. Their total number is 
fL/2i)' "where [x] is an integer part of a;. Examples of Ballot paths are shown on pictures a), 
b) and c) on Fig.l on page 7. Paths shown on pictures d) and e) are not the Ballot paths. 

Dyck paths are usually defined for L even and they are fixed at both ends as ho = Hl = 0. 
For L = 2p one has Cp := -Ti[ ^) Dyck paths which is the p-th Catalan number. For L odd 
close relatives of Dyck paths are those whose endpoints are fixed as ho = 1, h^ = 0. These 
paths are in one to one correspondence with the Dyck paths of length L + 1 and later on we 
will also refer them as Dyck paths. Among the paths shown on Fig. 2 on page 8 cases a), 
b) and c) are the Dyck paths, while cases d), e), f) are not. 

3 Raise and Peel Models with different boundary terms. 

3.1 The models definition. 

First, we describe the models algebraically and then, we discuss their physical interpretation. 

By definition, the ideal Xl forms left representation space of the boundary extended TL 
algebra. Consider on this space a dynamical process 

J^lPLit)) = -HLlpLit)), (3.1) 

Hl := E(l-e.) + c(l-/o)+c(l-/L), (3.2) 



defining an evolution of element \pL{t)) G X^. Here H^, the Hamiltonian of the process 
contains two numeric parameters — c and c, while the process itself depends also on four 
parameters of the algebra — q, a, a and b (see Eqs. (2.1) - (2.5) and (2.7)). We are 
interested in case where the boundary extended TL algebra becomes semigroup (that is, 
all the nonvanishing structure constants of the algebra are units) and so we fix algebra 
parameters as 

q = exp{i7r/3) (^ q + q~^ = 1), a = a = b = 1 . (3.3) 

In this case the Hamiltonian (3.2) becomes an intensity matrix and the process (3.1) can 
be given a stochastic interpretation (see, e.g., [10]). Expanding element \pL{t)) into linear 



combination of words of the ideal 

\PL{t)) = J2 PL{w\t)w 

weiL 

one treats coefficients Pl{u]\'^) ^^ unnormalized probabilities to find the stochastic system in 
configuration w at time t. 

In this paper we consider stochastic processes (3.1) corresponding to three particular 
choises of parameters c and c of the Hamiltonian (3.2). We call them models A, B, and C, 

model A: c = c = ; model B: c = 1 , c = ; model C: c = c = 1 . (3.4) 

In cases A and B the Hamiltonian acts invariantly on the subspaces of II spanned, respec- 
tively, by all Dyck and Ballot paths. Therefore we shall treat models A/B on their respective 
irreducible spaces of Dyck/Ballot paths. 

Now let us discuss physical interpretation of the models. We consider three processes of 
growth of a film of tiles which are deposited on a one-dimensional substrate of size L. 

As a substrate in all cases we choose profiles which are shown hatched on pictures (2.9). 
A rarefied gas above the substrate contains tiles and (possibly) half-tiles. They are moving 
along integer vertical lines as illustrated on picture (2.8) and upon hitting the substrate they 
can be absorbed and form interface configurations as shown on picture (2.9). Depending on 
a composition of the gas one distinguishes three cases 

model A: the gas contains tiles moving along lines with coordinates i = 1, 2, . . . , L — 1; 
possible interface configurations are given by Dyck paths; 

model B: the gas contains all the tiles and the half-tile moving along 0-th line; possible 
interface configurations are the Ballot paths; 

model C: the gas contains all the tiles and the half-tiles on both left and right bound- 
aries; possible interface configurations are the Anchored Cross paths. 

To determine evolution rules in the models we use the graphical presentation of the 
boundary extended TL algebra. We remind that in this presentation the substrate of size L 
corresponds to the unnormalized projector Xl; interface configurations correspond to words 
in the ideal Xl; the (half-)tile on i-th vertical line is an equivalent of the algebra generator 
Cj (/o//l for ^ = 0/L); hitting the interface by (half-)tiles amounts to left multiplication by 
Ci ifo, fh) in the ideal. With these identifications equation (3.1) defines following evolution 
rules. 

During an infinitesimal time interval dt a single event may happen with equal probability 
rate at any integer point of the interface. The following events are possible. 

a) At a local minimum point i (that is, if hi < hi±i) the interface either absorbs (half-)tile 
{hi I— > /ii+2) with probability dt or it refiects (half-)tile {hi stays unchanged) with 
probability 1 — dt. For the model C there is an exception from this rule described in 
item d). 
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a) 



b) 



c) 



d) 




Figure 1: The interface profile before event and the (half-)tilcs hitting the interface are drawn in dashed 
hnes. The interface profile after the event is drown in permanent line. Pictures a) and b) illustrate absorption, 
respectively, in a bulk and at the boundary of the interface. Pictures c) and d) are examples of avalanches, 
respectively, in the bulk (number of desorbed tiles n^ ~ 3, size of a substrate L = 6) and near the boundary 
(rid = 5, L = 6). Picture e) shows the total avalanche {nd — L — b). 



b) At a local maximum point i (that is, if hi > hi±i] 
and stays unchanged. 



the interface always reflects (half-)tiles 



At a bulk slope point i (that is, < ? < L and either /ij_i < hi < /ij+i, or hi-i > hi > 
hi+i) dropping a tile leads with probability dt to a nonlocal desorption event called 
avalanche. To describe the avalanche one determines integer k such that for all integers 
j standing between i and k inequality hj > hi holds and either hk = hi, or k runs out 
the interval [0,L], i.e., k equals L + 1, or —1. The avalanche causes desorption of one 
tile at each point j between i and k, that is hj t— *> hj —2. The avalanche size (a number 
of the desorbed (half-)tiles) nd = \i — k\ — 1, 1 < nd < L — 1, measures non-locality of 
the event. 

With probability 1 — dt the tile is reflected and the interface stays unchanged. 



= L that is. 



d) In the model C at a global minimum point i such that hi = 1 and hj > h 
dropping (half-)tile with probability dt causes total avalanche of a size nd 
hj \-^ hj — 2 , \f j ^ i. Again, with probability 1 — dt the (half-)tile is reflected and the 
interface stays unchanged. 



Typical absorption and desorption events are illustrated on Figure 1. 

For the model A the evolution rules described here were formulated in Ref. [12]. This 
stochastic process was named raise and peel model (RPM) there. Models B and C are 
versions of the RPM supplied with additional boundary terms. 



3.2 Stationary states: largest and smallest components and nor- 
malization factors. 

From now on we will study stationary states of the stochastic processes (3.1)-(3.4), i.e. 
solutions of equation 

Hl\pl) = 0. (3.5) 

Note that the intensity property of matrix Hi guarantees an existence of at least one non- 
trivial solution of equation (3.5). 








b) 



c) 



d) 



f) 



Figure 2: For the model A the minimal possible coefficient 1 in null eigenvector stands for pyramid type 
configuration a) for L even and for configurations b) and c) for L odd. For the model B half-pyramid 
configurations d) has coefficient 1 both for L even and L odd. For the model C the minimal coefficient nv^ 
appears again for configurations d) and for c)/f) in case of L even/odd. Dashed lines show the substrate in 
each case. 



Below we present results of a numeric investigation of the stationary states of three 
versions of RPM. Calculations were carried out with the use of REDUCE program for the 
system's size up to L = 13/11/10 for the models A/B/C, respectively. In all three cases the 
RPM has a unique stationary state. 

Denote components of \pl) (3.5) in the models A, B and C as p^ (w), Pi (w) and Pl {w), 
respectively. Here argument w labels in each case relevant interface configurations: those are 
the sets of Dyck paths {wj^j/cfc in the model A, Ballot paths {w} Ballot in the model B and 
Anchored Cross paths {w}ACross in the model C. Due to intensity property of the Hamilto- 
nian one always can choose null eigenvectors \pl) in such a way that all their components are 
nonnegative real (see [22]), thus, making the probabilistic interpretation consistent. More- 
over, it turns out that no one of coefficients Pi {w) vanishes. So, we can normalize them to 
be mutually primitive positive integers. Denote their smallest and largest components as 



m^L '■= niin{p5;'^(tf;)}. 



M 



■= raa,x{pl (w)} . 



(3.6) 



(a) 



m 



(b) 



1, but ml 7^ 1. The corresponding interface configurations 



It turns out that m^ 
are shown on Fig. (2). 

Denote total normalization factors of the stationary probability distributions as 



S 



(*) ._ 






(3.7) 



(*) (c) 

In the table below we collect values of St and m} for L < 11 



L 


1 2 3 


4 


5 


6 


7 


8 


9 


10 11 


4"' 


1 1 2 


3 


11 


26 


170 


646 


7429 


45885 920460 


4^) 


1 2 6 


33 


286 


4420 


109820 


4799134 


340879665 


42235307100 8564558139000 


4=) 


2 6 66 


858 


48620 


1427660 


47991340 


11589908610 


13642004193300 


1139086232487000 


-i^' 


1 1 2 


3 


11 


13 


10 


34 


323 


133 



(3.8) 



The table illustrates clearly the following three conjectures. These conjectures were found 
also in [2, 17]. 

Conjecture 1. Let A^ (resp., A^^ ) denote a number of vertically symmetric (resp., 
vertically and horizontally symmetric) alternating sign matrices of a size n x n ^. Then 

4? = ^r,« := (-3)"' n 2'^+ili - 4" = ^5S,. (3.9) 

l<i<p 

i<i<2p+i 

Conjecture 2. Looking at numbers standing in down-up diagonals of the table one observes 
equalities 

^(.) ^ si^^Sii,, (3.10) 

^P = rn^^S^l. (3.11) 

The first equality (3.10) allows one to get two expressions for S'2p_i ^ 

Ma) _ ,VH j,V _ ,VH ,,V ._ TT (3^ + 1) (6^) 1(2^)! 

•^2^-1 — ^4p-l/^2p-l — ^4p+l/^2p+l •— 11 {Ai)\{Ai + l)\ ' y^-^^! 

Their consistency, in turn, is based on relations conjectured in [23] 



aVH iaVH _ aV I aV _ V f l\2v-\j /o i o\ 

^4p+l/^4p-l ~ ^2p+l/^2p-l ~ , -|x/4p-2\ ■ \O.LO) 



(4p - 1) 



2p-l 



The second equality (3.11) relates quantities Sl and m^ . Looking at the table (3.8) one 
can also assume that vfvl are divisors of S^ . A formula for wtiJ refining this observation is 
guessed in [17] (see Eqs.(32), (33) there) 

my = Numerator of , s . (3.14) 



^On enumeration of various symmetry classes of alternating sign matrices see [23, 24] 
^Note that S2p_i is just the number of cyclically symmetric transpose complement plane partitions in a 
(2p)^ box (see, e.g., [25], p. 199). It is usually denoted as Ns{2p). 




N/\/\/\ [ \/\/ \^^ 

k, k=l,3,5,...,L. L 

a) b) c) d) 

Figure 3: In the model A the maximal coefficient in null eigenvector stands for configuration a) for L even 
and for configurations d) for L odd. In the model B configurations a) and b) enter the null eigenvector with 
the maximal coefficient if L is even and configuration c) corresponds to maximal coefficient if L is odd. 

Remark. At this point one may propose another reasonable normalization for the stationary state in 
the model C. Denoting components of the stationary state in this new normalization as {p^} one may fix 
normalization demanding that 

m^"^= min {j^^'> (w)} ^ si^\ (3.15) 

Certainly, the initial normalization is a more economic one and thus it is better suited for calculations. How- 
ever, normalization (3.15) has an advantage in interpreting results. E.g., formula for the total normalization 
factors in this case looks as 

q(c) _ q(a.) cib) _ 0(6) c(") _ 0(0)0(1) c(") /"q IK^ 

which is much more in the spirit of eq.(3.10) then the formula (3.11). We decide to keep initial 'economic' 
normalization throughout the text and to comment on the second normalization when presenting results. 

The last conjecture of this section describes the largest components of the stationary 
states in the models A and B. Analogous results for the model C are given in conjecture 13. 

Conjecture 3. The maximal coefficient MJf appear in the set {Pl } with multiplicity 1 for 
L even and (L — l)/2 for L odd. The maximal coefficient M| appear in the set {p^} with 
multiplicity 2 for L even and 1 for L odd. Their corresponding interface configurations are 
shown on Fig. (3) and values of the maximal coefficients are given by formulas 

M^,, = {S^'f , Mt.,^, = S^^l.S^l . (3.18) 



3.3 Model A: detailed distributions and Pascal's hexagon relation 

Let us first introduce more notation. 

An integer point i of an interface configuration {hk}k=o,...,L such that < i < L for the 
model A; < i < L for the model B; or, < ? < L for the model C, is called an N -contact 
ii hi = N and hi±i = A^ + 1. The contact points are local minima of the interface at which 
it can absorb (half-)tiles. For instance, the interface shown on picture d) of Fig. (3) has only 
one 1-contact i = 2 in case A, but it has two 1-contacts i = 0, 2 in cases B, or C. In all cases 
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it has 0-contacts at points k,k + 2, . . . L — 2 and additionally it has 0-contact at point L in 
case C. 

For either of the families of Dyck, Ballot, or Anchored Cross paths we denote {w}'^ 
(here * stands for an appropriate family name) the subsets of all interface configurations 
which have no p-contacts with p < N. These are the configurations whose global minimum 
(excluding points and L for the model A, or point L for the model B) is higher or equal 
to A^. Obviously, one has {wj^ = {w}^ and {w}^'^^ C {w}'^ . We shall stress that there is 
no strict correlation between the subsets {w}'^ corresponding to different sets of paths. For 
example, for any A^ 



i^JDyck <^ {^iBallot 5 {''^ I Dyck ^ {"^JBallot — 



where /c = 0/1 for L even/odd. This is because the left boundary point ho = 0, or 1 of a 
Dyck path is treated as 0-, or 1-contact point in a family of Ballot paths. 

Considering the subsets {w}^ proves to be useful for analysis of the raise and peel model 
in cases A and B. Here label A^ spans integers from to [{L — l)/2] in case A and from to 
(L — 1) in case B.'^ The subset with highest possible A^ in both models A and B contains only 
element which is drawn on Fig. (2), pictures a)/b) and d), respectively. As it was mentioned 
before, this element enters the stationary state with a minimal coefficient 1. In the sequel 
we determine detailed distributions S}^ j^ and detailed maxima M| j^ for each of the subsets 



N 



s 



(*) 

L,N 






(*), 



w] 



M, 



(*) 



L,N ■- 



m.ax{pi^' (w)} . 



(3.19) 



In this section we concentrate on studying the model A. The table below contains values 
of Si'j^ for L up to 13. 



L\N 


-1 





1 


2 


3 


4 


5 


6 


1 




1 
















2 


1 


1 
















3 




2 


1 














4 


3 


3 


1 














5 




11 


4 


1 












6 


26 


26 


5 


1 












7 




170 


50 


6 


1 










8 
9 


646 


646 
7429 


QK 


J 


1 

~8 










oo 

1862 


133 


1 








10 
11 


45885 


45885 
920460 


Af^ns - 


lOfi 


9 
276 




1 






4:OUO 

202860 


9660 


10 


1 




12 


93046J 


9304650 720360 18900 


375 




11 


1 




13 




323801820 64080720 2184570 


34452 


495 




12 


1 



(3.20) 



^By contrast, for the model C label N can take only two values — and 1. 
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For later convenience we add to the table a column for A^ = — 1 setting S'}z2p,Af=-i •= 

Sl=2pN=q- We also align data corresponding to even, or odd values of L leftwards, or 

rightwards in the columns, correspondingly. The numbers S')^"]y=o = ^l were already listed 
in table (3.8). 

It is helpful to look at the up-down diagonals in this table. The rightmost diagonal 
contains units only, the next one contains numbers (L — 1). For numbers staying in diagonals 
from third to seventh one finds expressions 



{L-2){L-3)-(2L+l) 

1 
2-3 

{L-2){L-3){L-4)(L-5)-(2L+l)(2L+3) 

22-32-5 

{L-3){L-4)2(L-5)(L-6)(L-7)-(2L-l)(2L+l)(2L+3)(2L+5) 

24.33.52.7 ' 

(L-3)(L-4)(L-5)2(L-6)2(L~7){L-8){L-9).(2L^l)(2L+l)2(2Z,+3)(2L+5)(2L+7) 



20.31.5^.7^.9 



26.34.53.72, 

(L-4)(L-5)2(L-6)2(L-7)2{L-8)2(L-9)(L-10)(L-ll).(2L-3)(2L-l)(2L+l)2(2L+3)2(2L+5)(2L+7)(2L+9) 

29. 35.54. 73. 92-11 

With these data one can guess general formula for Sf^]^. 



(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 



Conjecture 4 (Model A: detailed distributions). Denote n := [^^] — N . Integer n 
starts from and labels leftwards the up-down diagonals in the table (3.20). One has 



\ n — 1 1 r n — '2 i 



oia) _ 2^1.2/, " 1 ' A ' iL-[^]-p-l)\ Yr' (2L + 2n-6p-3)!! 

M (2P - 1)" L\ iL-2n + 3p)\ l[ (2L - 2[^] + 4p + 1)!! ^''''^ 

The same sequence was guessed in [17] ^. Our notation is related to that used in [17] as 
'^L!!v=[(L-i)/2]-n = R{n+l,L + 1) (see eq.(17) there). 

As it is explained, e.g., in [11, 17] the numbers given by formula (3.26) appear in counting 
elements of certain families of vertically symmetric alternating sign matrices. It is not 
however clear how one can characterize these families of ASMs. The table (3.20) suggests 
another combinatorial interpretation of integers Sj^j^. 

Pascal's hexagon relations. The numbers in sequence S^}^ (3.26) satisfy equalities 

c(n) c(a) I q{a) q{a) _ q{a) n{a) -r j ■ /o cy^\ 

'~'L~l,N'~'L+l,N+l^ '~'L,N-l'~'L,N+l ~ '-' L-l,N+l'~' L+l,N ■• ij 1^ iS euen , l-J-^'J 



^For a particular case L — 2p, N — 1 distributions 52^1 coincide with Pp{l), where Pp{k) is the unnor- 
mahzed probabihty to have k clusters in the stationary state (see [12]). An expression for Pp{k) is given in 
Conjecture 3 in [11]. 
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L=9 L=10 

a) b) 

Figure 4: In the model A configurations a)/b) of a size L = 9/10 have no 0- and 1-contacts. Thus, they 
belong to subset {wj^zfj.. In this subset among all the configurations of a fixed size L = 9/10 they have 
maximal possible number n = 2 of 2-contacts and in accordance with the statement of Conjecture 5 they 
enter the stationary state with a maximal coefficients -^fZo/in n-2' 

'-'l-i,n'-'l+i,n+i ^ '-'l,n'-'l,n+2 — '-'l-i,n+i'-'l+i,n ^ ij 1^ lb oaa. yo.zo) 

Substituting label N by n := ([^^^] — N) in the notation S^j^ one can uniformly write rela- 
tions (3.27), (3.28) as 

*-'L-l,n*-'L+l,n "r '-' L,n-1'-' L,n+1 ~ '-' L-l,n-l'-'L+l,n+l ■ {O.Z,U) 

In table (3.20) numbers participating in eq.(3.29) form a hexagonal structure. One such 
hexagon is drawn in the table. The components of each square monomial in (3.29) occupy 
opposite vertices of the hexagon. Relation (3.29) thus looks as a sophisticated variant of 
Pascal's triangle relation, wherefrom our notation follows. 

Like in Pascal's triangle case relations (3.29) can be used to reconstruct all numbers S*)^ „ 
provided some, say boundary, part of them are fixed. For certain boundary data the num- 
bers S'}"„ turn out to be integers and one recognizes among them the numbers of vertically 
symmetric and half turn symmetric alternating sign matrices, A^n+i ^"^^ ^n^ > whereas the 
total number of ASMs An, and the number of vertically and horizontally symmetric ASMs 
A^^ appear in combinations. A brief discussion of this subject is presented in the Appendix. 

We finish the section with observations of detailed maxima and of symmetry properties 
of the stationary distribution. 

Conjecture 5 (Model A: detailed maxima). In the stationary state of the model A 
maximal coefficient Mj^jy in the subset {w}'^ i^ stands for configuration with a maximal 
possible number {[{L — l)/2] — A^) of N- contacts (see explanations on Fig. (4)). It is given 
by formula 

^'^L,N = '^L-l,N~e{L) ' 

where e(L) := ~^^ ' is a parity function. ^ 

The relation (3.30) was also observed in [17]. Note that in case A^ = and L odd we use 
in eq.(3.30) an extension of sequence 5'jf]v ^o case A^ = — 1 that we made in table (3.20). 



<^ = 4"^i,^-ea) , (3.30) 



^Here and everywhere below in the main text we use definition of S^"}^ as it is given in eq.(3.26). This 
notation is natural for the stochastic models we are treating. Different convention about the second index 
in S*}^ „ is used in eq.(3.29) and in the Appendix. It is more suitable from a mathematical viewpoint. 
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Note also that equality Mi^J^2p+i,N=o = ^L=2p+i,N=i following from (3.30) agrees with the 
statements made in Conjecture 3 (see Fig. (3), picture d), cases k = 1 and k = L). 

Conjecture 6. In the model A for any configuration {hi} there exists a configuration {h'^} 
which appears in the stationary state with the same coefficient. The definition of {h'^} is 
following 

- for L even, h[ := /il-j/ 

- for L odd, let iq he coordinate of the leftmost 0-contact point in configuration {hi}, or 
Iq = L if there is no 0-contacts in {hi}. Then, h[ := hL_i + 1 for all i < L — iq, and 
h[ := hi-i — 1 for all i > L — Iq. 

For L even the statement of Conjecture 6 is a direct consequence of the left-right sym- 
metry of the Hamiltonian and the configuration space. For L odd the space of the states 
is no more left-right symmetric, and the symmetry is not an obvious one. Typical pair 
of symmetric configurations is drawn below (pieces of paths A/B and A'/B' are left-right 
symmetric). 



3.4 Model B: left orbits and factorizations to model A 

Again, we start with a few definitions. 

Consider a configuration w = {hk}k=o,...,L from the set of Anchored Cross paths. For any 
A^-contact point i of the configuration w we construct new configurations Wi{i) and Wr{i) 
(here value of A^ is irrelevant) 

r\ fu'\ w I hk + 2 for 0<k<i, , . 

Wi{i) = {h,}: h,:=\ , w./^^r (3-31) 



hk for i < k < L 



It may happen that either one or both configurations wi{i) and Wr{i) defined by eqs.(3.31), 
(3.32) do not satisfy condition c) of the definition (2.10) of Anchored Cross paths. Such 
configurations wi{i)/wr{i) should undergo total avalanche (see the end of Sec. 3.1), i.e., they 
are to be redefined as 

{h'J"}k=o,...,L -^ {h'J" - 2}fc=o,...,L , (3.33) 

Thus defined configurations wi{i) and Wr{i) are called, respectively, left and right coverings 
of the configuration w at the point i. Below, using the notion of left and right coverings we 
give an iterative description of left and right orbits of the configuration w. 

Let L{w, N)/R{w, N) denote a set of all A^-contacts of w which do not occupy right/left 
boundary point and which are placed to the left /right of any /c-contact with k < N. First, we 
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construct left/right coverings Wi{i)/wr{j) of w at all points i G L{w,N)/j G R{w,N). Note 
that for all coverings wi{i)/wr{j) the sets L{wi{i),N)/R{wr{j), N) lie strictly inside the set 
L{w, N)/R{w, N). Next, we construct left/right coverings for all obtained at a first step con- 
figurations Wi{i)/Wr{j) at all points from their corresponding sets L{wi{i),N)/R{wr{j),N). 
We repeat the procedure until at some step of iteration there will be no more A^-contacts 
available for generation of new coverings. The set of all thus obtained left/right coverings 
together with w itself is called left/right N-orbit of the configuration w and is denoted as 
Oi{w,N)/Or{w,N); w is called generating element of the orbit. The process of iterative 
construction of left and right orbits is illustrated on Fig.6, p. 30 and Fig. 7, p. 31, respectively. 
One can see that the numbers of elements in the orbits Oi{w, N)/Or{'w, N) and in the sets 
L{w,N)/R{w,N) are related as 

i^Oi{w,N) = 2#^("''^), i^Or{w,N) = 2#^('"'^). (3.34) 

An A^-orbit is called maximal if it is not a subset of some bigger A^-orbit. 

For the family of Ballot paths only the procedure of left covering makes sense and, hence, 
only the left orbits are defined. In the rest of this section we apply notion of the left orbits 
to establish detailed relations for the stationary states in case B raise and peel model. The 
maximal 0- and 1-orbits will be especially important for us. The reasons are following. 

The family of Ballot paths splits uniquely into collection of all (mutually nonintersecting) 
maximal left 0-/ 1-orbits. 

The maximal left 0-orbits of length L are generated by elements 

TA 1 u / for L even, 

w = {h,},=o,...,L : ^0 = I 1 fo^ ^ ^^^ (3.35) 

The maximal left 1-orbits of length L are generated by elements 

... . f either 0, or 2 for L even, /o o^\ 

w = {h,},=o,...,L : h, = l ^ for L odd. ^^'^^^ 



We further observe that sums of the stationary coefficients over maximal left 0- and 
1-orbits in case B are related to certain stationary coefficients in case A. These relations 
generalizing formula (3.10) are presented in 

Conjecture 7. Take any Dyck paths u and v of the lengths L and {L + 1), respectively. 
Denote *v a Ballot path of the length L which is a (left) reduction of the Dyck path v, i.e., 

hk(*v) = hk+i{v) , for all k = 0,...,L. 

With these notations one has 

fori even: ^ pf (^) = ^it p?(«) , E #M = ^i"^p£iW. (3-37) 
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fori odd: Y: Pfi^) = St^pf^M. E rfV) = Stlvfiu). (3.38) 

Here summation is taken over all elements of the left 0-, or 1-orhits generated by u, or *v 
(note that configurations u and *v span all generating elements described in (3.35), (3.36)). 

Combining the last conjecture with conjectures 4-6 from the previous section one can 
get exphcit values for sums over certain 0- and 1-orbits. These results and similar formulas 
for higher orbits are presented in the next 

Conjecture 8. Let us denote W{hQ., s) = {/iA;}A;=o,...,L (^ Ballot path of a size L such that 
for a given left boundary height Hq it has maximal possible number of s-contact points. We 
are interested in configurations where s take values between maxjO, {Hq — 1)} and (-^^i^ — 1) 
(see Fig. (5)). 

One can guess explicit expressions for sums of the coefficients p\j over left Hq- and 
{Hq — l)-orbits generated by elements W{ho,s). For even size L parameter ho takes on 
even values from to L and one has relations 

{ho = 2m) -orbits: ^ rfV) = S^Lli,mSt-i,s-mi 



v<^Oi(W (ho=2m,s),2m) 



(b)t,A _ cW c(a) ^^'^^^ 



{ho- I = 2m- l)-orbits: ^ p\^\u) = 5i"^_i ^l,,_^ 

ueOi{W{ho=2m,s),2m~l) 

The last formula is applicable in case m = 0, i.e., for (-l)-orbits. In this case all the 
(-l)-orbits are treated as singlets containing their generating elements only. 

For odd size L parameter ho takes on odd values from 1 to L and one has relations 



{ho-l = 2m)-orbits: J2 PlM = ^^L ^ti-™ , 

ueOi{W{ho=2m+l,s),2m) 

{ho = 2m + 1) -orbits: ^ pJ'V) = ^iti,™ '^Pi.^-m-i • 

ueOi{W{ho=2m+l,s),2m+l) 



This Conjecture is graphically illustrated on Fig. 6 on page 30. 

Most of the orbits appearing in the left hand sides of relations (3.39) and (3.40) are 
singlets and doublets. Namely, one has singlet {ho — l)-orbits in case s > ho and singlet 
/iQ-orbits ior s = ho — 1; one has doublet /iQ-orbits in case s > ho + 1. The larger size 
2~2— -plet ho- and {ho — l)-orbits appear, respectively, for s = ho and s = ho — 1. 

Relations for singlets and doublets provide explicit expressions for a number of stationary 
state coefficients in the model B. In particular, formulas for {ho — l)-singlets generated by 
elements W{ho, s) with s > ho were observed in [17]. 

Noticing that two configurations from /iQ-doublet generated by W{ho, ho + 1) correspond 
to {ho — 1)- and (/io+2)-singlets generated by W{ho, ho-\-l) and W{ho+2, ho + l), respectively, 
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W{ho,s) :-- 





b) 



c) 



Figure 5: A typical configuration W{ho,s) is shown on picture a). It is uniquely defined by its size L, 
its left boundary height Hq and the height s of its contact points in the bulk. A number of s-contacts 
in such configuration equals ( 1^'"" — s — 1). Boundary cases with s = {ho — 1) and s — { ^^^° — 1) are 
shown, respectively, on pictures b) and c). In the last case one uses parameter s just by analogy, because 
configuration c) contains number of contacts in the bulk. 

one concludes that expressions from the right hand sides of relations (3.39) and (3.40) should 
comply certain consistency rules. These are just Pascal's hexagon relations (3.27), (3.28) 
and in this way Pascal's hexagon relations were first time observed. 

We finish the section with several observations about detailed normalizations and ex- 
tremes in the model B. 

Conjecture 9 (Model B: detailed distributions I). For the quantities S]^)^ introduced 
in (3.19) one has formulas 



for L even: S\^ 



(b) 



N 



for L odd: S 



(b) 
L,N 



S 



S 



(a) 
(a) 



s 



(a) 



i-ffl 



s 



(a) 



(3.41) 



[^] 



The same formulas for 5']^jv were observed in [17] (see eq.(48) there). 

We remind that label A^ refers to global minimum of Ballot paths contributing to distri- 
bution S\^j^ (it should be greater or equal to A^). A^ may take all integer values from to 
L. For A^ = relations (3.41) coincide with eq.(3.10). 

One can consider another set of detailed distributions T,j^j^]^ which count stationary 
state contributions of all configurations w whose left boundary height is fixed as ho{w) = N 
and whose all s-contacts in the bulk satisfy condition s > M (i.e., whose global minimum in 
the bulk is greater or equal to M). For M > A^ all such configurations can be obtained by 
adding tiles (not half-tiles!) to the configuration W{hQ = N,s = M) drawn on Fig.5 a). 



Note that (in contrast to the case of S 



(b) 
L,Nj 



for distributions S^^^^^ label A^ may take 



only integer values between zero and L which are of the same parity as L. 



.(b) 



Conjecture 10 (Model B: detailed distributions II). In case M > N quantities T}j^'j^ j^ 
are given by formula 



<b) 



S 



(a) 
L+1,M~ 



[V 



Cf(a) 



(3.42) 



We were unable to find similar expressions in case M < N. 
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The right hand sides of eqs.(3.41) and (3.42) contain the same square combination 

) q(«) 



Si'K x^LY but for indices X, Y taking values in adjacent domains. Hence, formally speak- 



ing one can view distributions Sj;^j^_i and 5'}^jv-2 ^^ extrapolations of distribution S^-^vm 
to cases M = N — 1 and M = N — 2, respectively. However, we don't know any non-formal 
explanation to this fact. 

Conjecture 11 (Model B: detailed extremes). Among elements of the set {w}^^iig^ 
configuration with maximal possible number [^^^^] of N -contacts enters the stationary state 

with the maximal coefficient MJ^jy (^^^ def.(3.19)). All such configurations are generating 
elements W^(/io, s = Hq — I), or W^Hq, s = Hq), respectively, for Hq-, or {Hq — l)-singlet orbits. 
Hence, expressions for maxima M| jy '^'"c particular cases of formulas (3. 39), (3. 40) (index s 
there corresponds to N in the notation for detailed maxima). They are 

,. r ^^^L,N=2m — *-'L,m-l '~'L,m-> ,^ ,„< 

for L even: ,^^ , ^ , ^ (3.43) 

l\Ab) _ q{a) q{a) 

for L odd: ^^' . ^ . ^ 3.44 

In case of L even and N = there is one more configuration w = {ho = 2,hi = e{i)} (see 
Fig. 3 b) on p. 10) which enters the stationary state with the same maximal coefficient. 

Among all Ballot paths with the same left boundary height Hq configuration which have no 
any contact points in the bulk enters the stationary state with a minimal coefficient mih^- 
All such configurations are generating elements W^Hq, s = ■^+^" — l) of (Hq — I)- singlet orbits. 
So, expressions for m}^)^^ can be found among formulas eqs.(3.39), (3.40). They are 

^?M = 4j^] • (3.45) 



3.5 Model C: right orbits and factorizations to model B 

In this section we consider right 0- and 1-orbits in a family of Anchored Cross paths and 
describe corresponding stationary distributions for the model C. Our key observations are 
following. 

The family of Anchored Cross paths splits uniquely into collection of all (mutually non- 
intersecting) maximal right 0-/ 1-orbits. 

The maximal right 0-orbits of length L are generated by elements 

w = {hk}k=o,...,L ■ hL = 0. (3.46) 
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The maximal right 1-orbits of length L are generated by elements 
w = {hk}k=o ... L '■ either Hl = 0, or h^ = 2 and min hi = 1 (i.e., 7^ 0) . (3.47) 

' ' 0<J<-L 

We further observe that sums of the stationary coefficients over maximal right 0- and 
1-orbits in case C are related to certain stationary coefficients in case B. These relations 
generalizing formula (3.11) are given in 

Conjecture 12. Take any Ballot paths u and v of the lengths L and {L + 1), respectively. 
Denote v* an Anchored Cross path of the length L which is a (right) reduction of the Ballot 
path V, i.e., 

for all k = 0,...,L, hk{v ) = ■> 

[ hk{v) — 1 , otherwise. 



With these notations one has 

Oi \ CM/ X m-r Ota. 

qia) 



E #(z.) = Denominator of [-1^] p1\u) ^ ^^p^;:\u), (3.48) 



J2 Pl\i^) = Numerator of i—^\ p^Liiiv) = m^L^PLliiv). (3.49) 

ueOriV,!) \'Jl+2/ 

Here summation is taken over all elements of the right 0-, or 1-orbits generated by u, or v* 
(configurations u and v* span all generating elements described in (3.46), (3.47)). 

This Conjecture is graphically illustrated on Fig. 7 on page 31. 
Remark 1. In normalization (3.15) relations (3.48) and (3.49) read 

Y: pi^V) = 4tpf(-), (3.50) 

^ pi^V) - 4''^pSi(-)- (3.51) 

In this presentation Conjecture 12 looks very much in the spirit of Conjecture 7 and one can be immediately 
convinced that it refines formula (3.16). 

Remark 2. In the model C stationary state obeys mirror symmetry. Namely, any pair of configurations 
{hk} and {h'f.} whose shapes are left-right symmetric, i.e., 

h'l^ = /iL_fe for all k in case of L even; 

I ^L-fc + 1 , if miUj^Q hi — , 
h)^ — i T for all k in case of L odd; 

[ h-L-k - 1 , if miUj^o hi = 1 , 

enter the stationary state with the same coefficients. Using this symmetry one can reformulate Conjecture 12 
in the language of left orbits. 

19 



There are many singlets and doublets among the orbits appearing in the right hand 
sides of eqs.(3.48) and (3.49). E.g., all 0-orbits generated by elements w : hiiw) = 
and hk{w) >l\fk = l,...,L — 1, are doublets; all 1-orbits generated by configurations 
w : HlIw) = are singlets. Relations for singlets and doublets provide a number of explicit 
expressions for stationary state coefficients of the model C In particular, expressions which 
are presented in a right column of Table 2 in [17] all correspond to the right 1-singlets. 
To illustrate the practical use of Conjectures 12 and 8 we shall derive here formulas for 
stationary state coefficients of configurations 



X(s) := ho=2 




hr=2 



(3.52) 



even 



where configuration's size L may take even values only and parameter s runs from 1 to 
Denote 



{L-2) 
2 ■ 



Yis] 



ho=2 




Zis) := ho=3 



even 




According to Conjectures 8 and 12 the stationary state coefficients for these configurations 
and for configurations W{hQ = 0, s) and W{ho = 1, s) (see fig. 5 on p. 17) satisfy relations 

p'L\Yis))+p^^\x{s)) = !^^||p^ pP{Y{s)) , (s > 1) , 

'Jr. 



right 0-doublet (eq.(3.48)): 



right 1-singlet (eq.(3.49)): pt'^(r(s)) = m^ ptl,{Z{s + 1)) , 

left 0-doublct (first cq.(3.39)): P^l\W{0, s)) +P^l\Y{s)) = 4tl 4^-1,5 , («>!), 

left (-l)-singlct (second eq.(3.39)): Pl\W{0, s)) = ^["^ 5^ , 

left 1-doublet (second eq.(3.40)): p2'|i(Vr(l, S + 1)) + Pl+i(^(s + 1)) = S'^lU SQ , (s > 1) , 



left 0-singlet (first cq.(3.40)): 
wherefrom one obtains 



(a) r,{a) 



pZiimhs+1)) = 4t47w 



Jc) 



pf (X(5)) = m 



(c) 



' qia) qia) qia) ^ 

'^L+2 '^L+l '^L-l,s _ r, q(a) qia) qia) qia) I 

(a) ^ '^L+2 '-'L,s ^" '-'L+1 '^L+l.s+l ( 



(3.53) 
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In the last conjecture we describe largest components of the stationary states of the 
model C. 

Conjecture 13 (Model C: maxima). T/ie maximal coefficient Mj^ appear in the set {pi } 

with multiplicity 1 for L even and 2 for L odd. 

For L odd one of the corresponding interface configurations is the substrate shown on 
Figure 3 c) on page 10. The second one is a mirror image of the first (see Remark 2 to 
Conjecture 12). For L even the corresponding configuration is X{1) (see def.(3.52)). Explicit 
values of the maximal coefficients are 

for L odd: m]^^ = m^l^ S^^ Si%^, , (3.54) 



q{a) c,(a) c,(a) 

for L even: Ml' = ml' \ ^^ - SI Sll^^^ ) . (3.55) 



Relation (3.55) is a particular case of eq.(3.53) simplified with the use of Pascal's hexagon 
relations. 

A mirror image of the configuration X{1), which is the substrate, enters the stationary 
state with next to largest coefficient 

for L even: p I' {Substrate) = f^l Sl+i^l+ii- (3.56) 



4 Discussion 

Two of the issues discussed in this paper deserve further investigation. 

The first one is Pascal's hexagon relation which surprisingly comes as defining recurrent 
equation for some of the RPM's stationary coefficients. At the moment we do not have 
explanation to this fact. A very preliminary investigation of Pascal's hexagon relation is 
carried out in the Appendix. It shows that besides the solution S}^ p^ (3.26) which is related 
to RPM with open boundaries the relation admits also solutions which manifest themselves 
in the periodic stochastic models (see eq.(A.8) and the sentence above it; see also remark 
4 on page 24). Besides that, an extensive list of relations (A.27)-(A.38)) obtained in the 
Appendix indicates an existence of a profound relation between solutions of Pascal's hexagon 
reccurence and combinatorics of the alternating sign matrices. This relation has to be further 
explored. 

Secondly, while studying the stationary states of the models A, B and C we observed 
close relations between their probability distributions. Similar relations were also observed in 
Ref. [16] for the spectra of the models. Such relations have been interpreted by introducing 
the notion of orbits on the space of configurations of the models B and C. These facts 
suggest an idea that models B and A can be obtained by factorization of the model C 
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It would be interesting to find out such conjectural factorization at the level of boundary 
extended Temperley-Lieb algebras. 
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Appendix. Solutions of Pascal's hexagon relation and 
ASM numbers 

Consider a set of numbers which are placed on vertices of a planar trigonal lattice. 
Introducing coordinate lines on the lattice as it is shown on figure (A.l) we label the numbers 
in the set by a pair of integers: {fm,n}, m,n E Z. 



m (n — coTist) 

(A.1) 




n / (m^const) 

For any six numbers a, a', b, b', c, c' in the set which are situated on apices of a hexagon 

c b' 



a -V a' (A.2) 



b c 



we impose a condition 



aa' + bb' = cc' , (A.3) 

which can be viewed as a generalization of Pascal's triangle relation to the case of hexagon. 
In coordinatization (A.l) it reads 

Jm~l,n /m+l,n ~r /m,n— 1 Jm,n+1 /m— l,n— 1 /m+l,rt+l • \-^'^/ 
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As the author has learned from Yuri Stroganov Pascal's hexagon relations belong to a 
wide family of bilinear equations, including Somos sequences, number walls, S-arrays, cube 
reccurence etc., which are currently extensively investigated in algebraic combinatorics (see, 
e.g., [27] and references therein; see also bilinear forum on the web). Moreover, it has 
been mentioned in physical literature under the name of discrete Boussinesq equation (see 
eq.(8.11) in [26] and the references therein). In [26] it was obtained as certain 2-dimensional 
reduction of the Hirota's difference equation, which is also called the octahedron reccurence 
by combinatorialists. 

A particular solution of Pascal's hexagon relation is given by formulas 



[^1 



(m 



\ n — 2 1 

^ .„-., . ,... [ii±2]-p-l)! [:p:] (2m + 2n-6p-3)!! 

•^"'•" ^t (2P - 1)!! ;={ im-2n + 3p)! /j, (2m - 2[^] + 4p + 1)!! ' ^ ^ 



r Ti—'2 1 



^ .-[.^/4, ^1 1 'fr' (m-[^]-p-2)! ':p:' (2m + 2n-6p-5)!! 

'"" M(2P-1)" io {m-2n + 3p + 2y. IJ, (2m - 2[^] + 4p - 1)!! ' ^ ' 

where m is an integer and n is a nonnegative integer. Here for negative values of m ratios 
of factorials are treated as Pochhammer symbols: 

71-1 

(n - iy./{k - 1)! = n J =■■ {k\n-k) , (2n - l)\\/{2k - 1)!! = 2("-'=)(A;)(„_fc)(2A;)(2„_2fc) • 

j=k 

The corresponding distribution of numbers on the lattice is shown in table 1 on page 32. 

Formulas (A. 5) and (A. 6) can be reproduced in a following way. First, one assigns certain 
initial data to fm,n (like in Pascal's triangle case) 

Jm>o,n=0 = 1 ? Jm>0,n=-l = 1 5 Jm=2n-l,n>0 = , (A. 7) 

and calculates values of fm,n for n > and m > 2n row by row (i.e., consecutively for 
n = 1,2,3, etc.) using relations (A. 4). In this way formula (A. 5) for detailed distributions 
^Ln — fm=L,n was obtained in Sec. 3.3. 

Then, one observes that formula (A. 5) is well defined on a half-plane n > and gives 
numbers satisfying Pascal's hexagon relation in this sector. Note that combinatorial func- 
tions R{n,L) and Q{n,L) introduced in [17] (see eqs.(16) and (17) there) both are related 
to function (A. 5) 

/„,„ = i?(n+l,m+l), /_„,„= (-l)['^]g(n+l,m + n + l). (A.8) 

Note also that in the sector m < 2n — l, n < 2m — 1 numbers fm,n vanish (see tab.(l)) that, 
in particular, destroys the connection between numbers fm,n in sectors m > 2n, n > and 
2m < n, n > via Pascal's hexagon relations. 
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Further on, comparing numbers fm,n in sectors m >2n > and n > 2m > one observes 
TTi <-* n symmetry 

fm,n = (-l)(™"+ W'»+«)/3]) 2-[(™-2n)/3] f^^^ V m > n > . (A.9) 

Extrapolating this symmetry relation one defines numbers fm,n in sector m > > n through 
those from sector n > > m. The Pascal's hexagon relation for the newly defined numbers 
stays valid. 

Finally, one can guess general formula (A. 6) for /„ „ on a half-plane n < analyzing the 
numbers in sector m > > n. The numbers fm,n thus obtained satisfy relation (A. 4) on the 
whole lattice. Let us end up discussion of the solution (A. 5), (A. 6) with a list of remarks. 

1- fm,n take particularly simple values along the following directions (see tab.l) 

fmfl = fm,-l = 1 (A. 10) 

/o,„ = (-1)[^"/31 2["/3] , /o,_„ = (-l/2)[2"/3] , V n > , (A.ll) 

/_^_„ = (_2)[(n+2)/3] ^ f_^_^ ^ (_l)[(2-n)/3] ^^/2y(2n-l)/3] ^ Vu > 1 , (A.12) 

Jm=2n-l,n>0 = /m>0,n=2m-l = /m=2n+3.n<-2 = Jm<-2,n=2m+3 = U. (A.ioj 

Conditions (A.10)-(A.13) can be used as initial data for reconstruction of all the 
nonzero numbers fm.n by means of relation (A. 4). 

2. Symmetry relations (A.9) are valid for any pair m, n such that m > n (but not for 
m < n). 

3. The numbers fm,n given by eq.(A.5) are always integer. The numbers f-m-n given by 
eq.(A.6) become integer upon multiplication by 2*^""^^. This fact is a particular mani- 
festation of the Laurent property, which was proved in general in [27] (for specialization 
to the case of octahedron recurrence see [28]). ^° 

4. The numbers fm,n in sectors m > 2n > and < n < 1 — m are known to be related 
to open and periodic one dimensional stochastic models, respectively (see [17] and 
Sec. 3. 3). One may expect that there exists some one dimensional stochastic process 
related to the numbers fm,n in sector n < — 1, m < 2n + 2. 

Next, we consider solutions of Pascal's hexagon relation which are polynomials in one, or 
several variables. Evaluated at particular values of their variable (s) they produce numeric 
solutions to (A. 4). We found two examples of polynomial solutions. 

First one is a set of polynomials in two variables Fm,n{x,y) which are defined in the 
sector n > 0, m > 2n — 1. They satisfy Pascal's hexagon relation (A. 4) (just substitute 

-"^•^The author is grateful to David Speyer for pointing out this fact to his attention. 
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fm,n by Fm,n{x,y) there) and one can uniquely define them provided that polynomials on 
the boundaries of the sector are fixed. We choose 



m>0,0 



1, Fm>l,l = y + (m - 2)X , Fm=2n-l,n>0 = , (A. 14) 



which is a generalization of boundary conditions (A. 7). A nontrivial fact is that with such 
choice of the boundary conditions one obtains polynomial (rather then rational) expressions 
for Fm,n{x,y). In table 2 on page 33 some polynomials Fm^n{x,y) for small values of m and 
n are listed. 

The functions Fm^n{x,y) obey rescaling symmetries'^ 

FmA^^y) = l/"^m,n(a:/t/, 1), for y ^ , (A.15) 

F^,„(x,0) = a;"F„,„(l,0), (A.16) 

F2n,n{x,y) = yF2n,n-l{x,y). (A.17) 

It follows then that only two values y = and y = 1 are essential. We keep using variable y 
as it is suitable for producing integral solutions to Pascal's hexagon relations. In particular, 
one finds following relations between polynomials -Fm,n and the numbers given by formulas 
(A.5) and (A.6) 

-^m+l,?i(l, 0) = Fm+2,„+l(0, 1) = Fm,?i(l, 1) = fm,n, (A. 18) 

Another set of polynomials Gm,n{x) are solutions of a variant of Pascal's hexagon relations 

^m,n-i-l^Tn,n—l ~r LrjTt+l,?!— l^m— l,n+l ^m—l,n'~^m+l,n ■ yJ\.Z,\J) 

These relations result from a different coordinatization of Pascal's hexagon rules (A. 2), (A. 3). 

^^Rescaling properties (A.15), (A.16) are particular manifestations of "gauge" transformations on the set 
of solutions of Pascal's hexagon equation. Namely, once a solution {a,a',b,b',c,c'} of eq.(A.3) is given a 
family of gauge equivalent solutions {aa, a'a', j3b, /3'6', 7c, 7'c'} is parametrized by six scaling numbers a, a' , 
P, f3' , 7 and 7' satisfying constraints aa' = (3f3' = 77' 7^ 0. On a trigonal lattice this gauge transformation 
is uniquely defined provided one fixes 6 mutually independent scaling parameters for the vertices situated 
as shown on a figure 
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The coordinate frame for them looks as 




n (m — const) 

(A.21) 



m (n^ const) 

For this set of functions we impose boundary conditions 

Go,n>o = l, Gi,„>o = l + nx, G„>o,o = :r[('"+i)/3l(a: + l)['"/3](-l)['"/3]+[™/^ (A.22) 

which allow one to iteratively reconstruct polynomials Gm,n{x) (again, a nontrivial fact) in 
the sector m > 0, n > 0. In particular, one finds 

Gm,l = a;['"/3] (a; + l)[(™+2)/3] (_i)[(m+2)/3] + [(m+l)/2] _ ^j^^S) 

Boundary conditions (A.22), (A. 23) are generalizations of the boundary data for /„,« along 
directions m < 0, n E {0, 1}, and n > 0, m E {0, —1} (see table 1 and formulae (A.IO)- 
(A.12)). Explicit expressions for Gm,n(fl) for small values of m and n are given in table 2 on 
page 33. 

The functions Gm,n{x) obey reflection symmetry 

G^,n{x) = a;[('"-"+i)/3]^-[(m-n+i)/3]](_^)[(m+n)/2]^_(~)^ for x ^ 0, -1, (A.24) 
G'„,n(0) = (-l)'"+^G'„_i,„+i(-l), (A.25) 

where x = —1 — x, fh = n — 1 and n = m + 1. At certain values of x they reproduce 
numbers (A. 5) 

^m,nyJ ) Jm+n—l,m—l j *-^m,nl,-Lj (^ ij J—n,m- yA.ZOJ 

Finally, it is remarkable that the solutions of Pascal's hexagon relation reproduce in 
several ways numbers which were discovered in counting various symmetry classes of the 
alternating sign matrices (see [23, 24, 25]). Below we present list of such formulae. Rela- 
tions (A.27), (A.28) and (A.31)-(A.34) were observed in [17] (see eqs.(18)-(23) there). First 
equalities in (A.27) and (A. 29) are particular manifestations of eq.(A.17). 

-^2?i,n-l(l5 1) = -^2n,?i(l, 1) = G'„_|_i^„(0) = G'n,n+l(0) 

= {1, 3, 26, 646, . . .} = A^+i , n > 0, (A.27) 
i^2„-i,n-i(l,l) = G'„,„(0) = {1,1,2,11,170,...} = -f±i, n>0, (A.28) 

^2n±l 
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i^2n,n-i(2,3) = ^F2„,„(2, 3) = {1,7,143,8398,...} = 4^, ^ > 0, 

F2„+i,„(2,3) = {1,5,66,2431,252586,...} = ^-!^, n > 0, 

G„,„(l) = {1,2,10,140,5544,...} = AfJ, n > 0, 

G'„,„+i(l) = {1,3,25,588,39204,...} = Af„^,, n > 0, 

G'„+i,„(l) = {1,2,2x7,7x42,42x429,...} = A„A„+i, n > 0, 

G'„,n+2(1) = {l,2^7^42^429^...} = (A„+i)2, n>0, 

^"'"v 2 '- 2' 2^' 2^' 2"^' 2^ ' 2*^ " 2" ' "- ^ ^j 

1 _ 1^2x6 6x33 33x286 A^^U3^ „^n 

1 1 3"^ 26^ [0 forn > 1 odd, 

^«."+i(-2) = {1'0'^'0'^'0'^'---} = | (iLtil! forn>Oeven, 

1 _ 1 1^ l3x3 1x3^ 3^x26 3x26^ 

<-^ra+2,nl~2'' ~ 1-2' 22 ' 2^ ' 2^ ' 2^ ' 2^ ' ' ' 



A.29) 

A.30) 

A.31) 
A.32) 
A.33) 
A.34) 

A.35) 
A.36) 
A.37) 



= ^(AV2)/4]+l)"^"+'^4n/4]+3)^^"^ n>0, (A.38) 

where p{n) = 1, 3, 3, 1, if (n — 4[^]) = 0, 1, 2, 3, respectively. 

Here A„ is a total number of n x n alternating sign matrices. Notations A]^, A^^ and 
A^'^ stand, respectively, for numbers of vertically symmetric, vertically and horizontally 
symmetric and half turn symmetric n x n alternating sign matrices. Expressions for A}^n+i 
and A^^j^i/A^,^^^i are given in (3.9) and (3.12). Formulas for A„ and A^^ are 

_- M3A: + l)! .HT_., ^2^^3A: + 2 ^t _ A 4 A3A:)!A:!y 
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0-orbit, octet, sum over orbit = 170 x 11 

676 





676 



234 



130 




2-orbit, quartet, sum over orbit = 50 x 4 



4-orbit, doublet, sum over orbit = 6x1 





Figure 6: Examples of maximal left orbits for size L = 6 Ballot paths. 

Arrows on pictures point from source configurations to their coverings; bold dots indicate 
points where coverings happen; dashed lines show configuration before covering happen. 
Numbers standing below configurations are the corresponding stationary state coefficients of 
the model B. Sums of the coefficients over orbits are described by the first one of relations 
(3.39). 
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0-orbit, octet, sum over orbit = (m^ Sj /S^ ) x p^ 



(b) 



170 X 78 



6050 





4226 



625 



1430 




1-orbit, quartet, sum over orbit = m^ x pg 




11 X 676 




2286 



159 



Figure 7: Examples of maximal right orbits for size L = 5 Anchored Cross paths. 

Notations are the same as on Fig.6. Numbers standing below configurations are the corre- 
sponding stationary state coefficients of the model C. On the last picture two configurations 
suffer from the total avalanche (3.33). 
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m=2n+3 n=2m+3 





429 143 _22 

32 32 32 



/ / 



1 J_ 14 

8 16 32 



252586 26194 



16 16 

22610 8398 2431 429 

8 8 8 8 







66 18 1 1 12 84 







646 425 260 143 66 21_ 

4 4 4 4 4 4 



16 16 



'il 1 _i i 

8 8 4 4 



16 16 



5 49 588 

8 8 8 







_5 _1 1 _4 _25 _70 147 

4 2 4 4 4 4 4 

17151311.9 7 5 3 1 _1 _3 _5 _7 _9 _11 

2 22222222222222 



1 1 1 1 1 



1 1 1 1 1 1 1 1 

/ / 



1 1 1 1 1 1 1 1 1 1 



-8 -7 -6 -5 -4 -3 -2-1 1 



1 1 1 1 

2 3 4 5 



-195 -132 -84 -49 -25 -10 -2 11 

/ / 

4356 1764 588 140 14-2 2 2 



3 11 26 50 
^ 26 170 







39204 5544 294 -18 4-2 3 3 

/ / 
-18018 528 -44 4 2-8 -11 



-60 -8 -4 -12 100 170 



n=-l 



n=0 



, m=2n— 1 



-1716 52 4 4 -10 26 26 
/ / 



m= — 1 



m=0 



n=2m— 1 



Table 1: Particular solution (A. 5), (A. 6) of Pascal's hexagon relation (A. 4). The numbers standing 
along indicated lines are given by formulas (A.10)"(A.13). Knowing them one can reconstruct data 
on the whole lattice by a repeated use of Pascal's hexagon relation. 
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F4,2 = 2xy + y'^ ^5,3 = llx'^y + I2xy'^ + ?,y^ 

F5,2 = 3x2 + 6xt/ + 2|/2 ^7,3 = 26a;3 + 78a;2|/ + 55x1/2 + Iiy3 

F6,2 = lla;2 + 12a;i/ + 3t/2 ^§,3 = llQx^ + 2Mx'^y + l^Qxy"^ + 2Qy^ 

F7,2 = 50a;2 + SOajy + 5t/2 ^9,3 = 646a;3 + SIGx^t/ + 350a;|/2 + 50|/3 

^8,2 = 85a;2 + 42x|/ + Qy"^ 

F9,2 = 133x2 + 56a;i/ + 7y2 

F8,4 = 170x^1/ + 294a; V + 156x1/^ + 26|/^ 

F9,4 = 646a;^ + 2584x^1/ + 2839a;2i/2 + ngOxy^ + 170i/^ 



G2,2 = 2 + 5a; + 3x2 q^^^ ^ 3 + 73, + 43,2 

G2,3 = 3 + llx + llx2 ^3,3 = 11 + 44x + 59x2 + 26x3 

G2,4 = 4+19x + 26x2 ^3,4 = 26 + 137x + 255x2 + 170x 

G2,5 = 5 + 29x + 50x2 G'3_5 = 50 + 321x + 747x2 + 646x3 

Gofi = 6 + 41x + 85x 



.3 



.2 



^4,2 = -4X - 9x2 _ g^3 

G4,3 = 26 + 97x + 121x2 + 50x3 

^4,4 = 170 + 935x + 1956x2 + 1837x3 + 646x^ 

G5,2 = 5x +16x2 + 17x3 + 6x^ 

^5,3 = -50x - 179x2 - 214x3 _ §5^4 

^6 2 = -6x - 19x2 - 20x3 - 7x^ 

Table 2: Lattice polynomials Fm,n{x-,y) and Gm,n{x) for small values of m and n. 
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